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The evolution of iron catalyst particles during aerosol (floating catalyst) Chemical Vapor Deposi¬ 
tion (CVD) synthesis of Carbon Nanotubes (CNTs) from CO is computed using a multi species 
Computational Fluid Dynamics (CFD) model incorporating a lognormal aerosol method of moments 
(MOM) to describe their dynamics and a combined chemical kinetics and equilibrium model for cat¬ 
alytic production of C0 2 . The influence of the presence of iron at the reactor walls, the fed particle 
size, number concentration and polydispersity and the effect of the catalytic production of C0 2 at 
the reactor wall are studied in terms of particle size, concentration and polydispersity and reagent 
concentration during CNT synthesis. It is found that iron catalyst particle dynamics are essentially 
insensitive to wall iron concentrations and, for a wide range of particle sizes and concentrations, it 
is found that the catalyst particles are stable up to a critical CNT window in which CNT nucleation 
and growth occurs. Concentrations of catalyst particles significantly above 1 x 10 14 #/m 3 , however, 
lead to poor control over catalyst particle size and polydispersity at the CNT nucleation front which, 
in turn, leads to poor control over CNT diameter. The location of the growth window is, shown to 
be directly associated with the availability of catalytically produced C0 2 diffusing from the reactor 
walls to the reactor core. These results help to explain the large variations in CNT diameter and 
chirality and the inefficient use of catalyst material in other floating catalyst CNT processes based 
on in-situ catalyst particle synthesis. 

Keywords: Carbon Nanotubes, CNT, Computational Fluid Dynamics, CFD, Catalyst, Aerosol, 
Particle, Method of Moments, Evaporation, Coagulation, Population Balance, Chem¬ 
ical Kinetics. 


1. INTRODUCTION 

Carbon nanotubes (CNTs) are of interest due to their 
unique electrical, structural, chemical, thermal and optical 
properties, however, many applications of CNTs require 
large volumes of functionally similar material having a 
narrow range of properties such as diameter, length and 
chirality. For instance, Single-Walled CNTs (SWCNTs) 
can range from wide band gap semiconductors to metal¬ 
lic depending upon their chiral index . 1 It is known that, 
for aerosol (floating catalyst) Chemical Vapor Deposition 
(CVD) synthesis methods, the size of catalyst particles is 
critical in determining CNT diameter and thus, to some 
extent, the chirality and that residence times under growth 
conditions help to determine CNT length . 2 Consequently, 
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controlling catalyst particle sizes and residence times in 
growth environments are important requirements for ulti¬ 
mately controlling the properties of produced CNTs. Poor 
control of these parameters can lead to low efficiency 
in the use of catalyst material, production of significant 
amounts of carbon materials other than CNTs and wide 
distributions of CNTs properties. 

Traditionally, floating catalyst methods for CNT pro¬ 
duction have relied on in-situ formed catalyst particles 
produced from chemical precursors in the region of 
CNT nucleation and growth . 3,4 These methods produce 
CNTs having wide distributions of diameter, length and 
chirality and result in product containing a significant 
amount of unused catalyst material and amorphous car¬ 
bon which requires further purification steps. Recently, we 
have developed a novel aerosol method for the synthesis 
of CNTs relying on the introduction of well controlled 
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pre-made catalyst particles into an aerosol reactor via 
a hot wire generator (HWG). 5 It was shown that the 
method is useful for the production of both single-walled 5 
and multiwalled 6 CNTs with controllable dimensions of 
catalyst particles and CNTs. 2 In this work, we aim to 
understand the role of catalyst particle evolution in CNT 
synthesis by modeling the flow, heat transfer, gas phase 
species production and transport as well as particle trans¬ 
port and dynamics in our reactor. 

Other work on the reactor level modeling of CNT pro¬ 
duction can be found in the literature, however, no previ¬ 
ous models have yet been applied to predict polydisperse 
CNT catalyst particle dynamics and transport. The sim¬ 
plest reported modeling studies neglect the catalyst particle 
dynamics and model only the reactor flow and temperature 
conditions (e.g.. Ref. [7]). Scott et al. 8 have studied the 
HiPCO process 4 by using the results of Eulerian calcula¬ 
tions of the flow inside the reactor as temperature bound¬ 
ary conditions for trajectory calculations of iron cluster 
dynamics. Cluster calculations included detailed reaction 
mechanisms for a family of iron carbonyls and the growth 
of Fe„ and Fe„CO clusters for n up to 999. Importantly, 
Scott et al. 8 did not include cluster collision processes in 
their model. Endo et al. 9 and Kawana et al. 10 developed a 
CFD model that focuses on the chemistry of nanotube syn¬ 
thesis but did not consider the catalyst particle transport or 
dynamics in the reactor. 

The computation of catalyst particle dynamics in CNT 
reactors can be carried out using a variety of numerical 
techniques, but such phenomena are usually modelled by 
solving the flow and gas phase partitioning in an Eulerian 
reference frame and tracking particles through the flow 
field in a Lagrangian reference frame. For highly diffu¬ 
sive particles, statistical accuracy requires a large num¬ 
ber of trajectories and small time steps to account for 
Brownian motion. Significant reductions in computational 
effort can be achieved using an Eulerian formulation for 
multi-phase processes. In general, these can be classi¬ 
fied as monodisperse methods, discrete population bal¬ 
ance (DPB) (or sectional methods) and moment methods. 
Recently, Kuwana and Saito extended their CNT reactor 
model to include nucleation and growth of catalyst parti¬ 
cles assuming monodispersity via a two moment method. 11 
In this work, we use a three moment method in which the 
problem of mathematical closure is addressed by assuming 
the distribution is lognormal, 12 thus allowing the polydis- 
persity of the Particle Size Distribution (PSD) also to be 
computed. 

In this paper this validated model is used as a pre¬ 
dictive and diagnostic tool to understanding factors that 
influence the control of CNT synthesis in floating cata¬ 
lyst methods. In particular, we investigate the behavior of 
pre-made catalyst particles up to the location of the exper¬ 
imentally observed CNT nucleation front. The influence 
of iron vapor concentrations at the reactor wall and fed 
particle size, number concentration and polydispersity are 


computationally investigated and correlated with CO wall 
reaction kinetics and the transport of produced C0 2 . 

2. MODELED EXPERIMENTAL SETUP 

Our method for the production of SWCNTs has been 
described elsewhere. 5,6 In the method, pre-made catalyst 
particles of controlled size, concentration and polydisper¬ 
sity are introduced into well defined synthesis conditions. 
Briefly, in this method, Fe particles were produced from 
Fe vaporized from a resistively heated catalyst wire in 
a H 2 /N 2 (with a 7/93 mol ratio) flow (400 cm 3 /min). 
Here, particles were formed and grown by vapor nucle¬ 
ation, condensation and particle collision processes. Subse¬ 
quently, the produced particles were introduced into either 
a ceramic or stainless steel tubular reactor, mixed with a 
carbon monoxide flow (400 cm 3 /min) and heated to induce 
CNT formation. Reactor wall peak temperatures ranged 
from 1163 K to 1480 K. In our first experiments, iron 
and carbon monoxide were used as the catalyst material 
and the carbon source, respectively. Experimental obser¬ 
vations of CNT synthesis behavior in the reactor at a peak 
wall temperature of 1197 K concluded that there is a spe¬ 
cific “growth window” for CNT nucleation and growth in 
which all CNT synthesis occurred, however, CNT produc¬ 
tion was unstable. 13 Further analysis of the experimental 
data showed that these instabilities were primarily associ¬ 
ated with the reactor wall conditions. It was found that, in 
order to provide reproducible CNT synthesis conditions, 
the reactor walls had to be saturated with iron. This was 
achieved either by deposition of iron particles on the reac¬ 
tor walls or by using reactor tubes made principally of 
iron, i.e., stainless steel. Depending on the level of iron 
saturation at the wall, it was observed that CNTs did or did 
not form. Reactors with iron saturated walls (either via the 
stainless steel reactor wall tube or by preconditioning with 
iron vapor and/or particles) were found to produce CNTs 
while reactors with “clean” walls (unconditioned ceramic) 
produced no CNTs. At these conditions, the growth win¬ 
dow was determined to be between 1167 and 1179 K in 
the heating region of the reactor 13 when the wall reaction 
temperature was 1480 K. This temperature coincides with 
the results reported in Ref. [5], where the growth of CNTs 
in a ceramic tube was observed above a maximum reactor 
temperature of 1163 K. Similar temperature windows for 
CNT growth were found with CNT production via thermal 
decomposition of ferrocene. 14 From this data, we devel¬ 
oped the hypothesis that the iron vapor pressure at the 
reactor wall were critical in stabilizing catalyst particles by 
reducing or eliminating catalyst particle evaporation before 
CNT nucleation. 

In the current modeling work, boundary conditions are 
chosen to simulate conditions wherein reactor walls were 
saturated and unsaturated with iron at the highest tempera¬ 
ture reactor experimentally investigated in Ref. [13] (with a 
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Outlet 



Power supply 


Fig. 1. Schematic of modeled experimental setup showing the wall tem¬ 
perature profile and the location of the CNT “growth window.” 


maximum reactor temperature of 1480 K) where, it would 
be expected, particle dynamics would be most pronounced. 
A schematic of the experimental apparatus together with 
measured reactor wall temperatures specifying the location 
of the growth window is shown in Figure 1. Based on these 
conditions, a series of CFD calculations was performed to 
investigate catalyst particle dynamics and gas phase species 
composition. 

3. MODEL EQUATIONS AND BOUNDARY 
CONDITIONS 

Our previous work has modeled polydisperse aerosol 
transport in two and three-dimensions as well as the effects 
of chemical reaction and gas phase species mixing and 
strong heat and mass transfer coupling on the evolution 
of the aerosol size distribution. 15-19 This method uses an 
Eulerian moment form of the general dynamic equation 
for aerosol behavior developed using a lognormal size 
distribution in conjunction with a Reduced Navier-Stokes 
(RNS) formulation 20 for bulk fluid modeling with Mach 
number based streamwise pressure flux splitting and axial 
viscous terms included as a deferred corrector. 

Three sets of governing equations are considered: those 
related to the transport of mass, momentum and energy in 
bulk fluid flow, those related to the conservation of indi¬ 
vidual species in the fluid and those related to the transport 
and dynamics of suspended particles within the fluid flow. 
For all cases studied, the local Reynolds number based on 
the inlet diameter was below 500 and the Grashof number 
was greater than unity indicating laminar flow under the 
influence of buoyancy. Details of the formulation, coupling 
and implementation can be found in Refs. [17-19]. A brief 
summary, including additions to the model needed for the 
present calculations, are given below according to the vari¬ 
able definitions described in the Nomenclature Section. 


The bulk fluid flow is described by the compress¬ 
ible Reduced Navier-Stokes 20 equations with a deferred- 
corrector (DC) 21 to recover the Full Navier-Stokes 
solution. In non dimensional form, the governing equations 
with source terms due to fluid/aerosol particle coupling 
and with multiple species in the gas phase including the 
effects of buoyancy in each coordinate direction (pg_, 
can be written in vector form as 
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As, for the cases computed in the work, the mass fraction 
of particles is low, A H t is assumed to be zero. Closure of 
the equations is achieved with the equation of state. 


yocMi^-P = p^H,-^(y 00 -l)Mi(u 2 + v 2 + w 2 ) 
The total enthalpy of the mixture is defined as 

M 2 + l> 2 + W 


( 2 ) 
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- L - f - + h (3) 


according to the formulation given in. 18 

The present method uses an Eulerian moment form of 
the general dynamic equation (GDE) developed assuming 
a lognormal size distribution, 12,22 which has been found 
to exist in a wide range of aerosol systems. 23 The model 
incorporates contributions for particle transport (due to 
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convection, diffusion and thermophoresis) and dynamics 
(due to coagulation, evaporation and condensation). Nucle- 
ation is not considered in this work since the high con¬ 
centrations of catalyst particles introduced into the reactor 
provides sufficient surface area for condensation alone to 
maintain low supersaturation ratios throughout the reactor, 
thus eliminating the formation of additional particles. 16 In 
non-dimensional form, the GDE yields a moment form 
representing three simultaneous partial differential equa¬ 
tions, thus 


d(M k /p) 

p ^ +pU/' ■ V(M k /p) 


dt 


■ V- 
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Sc k Re t 


VM, +V- [M, 


K T 


; Re, T, 


- V7’, 


r.'f 


= y k M 2 k M koo t x + J2 ( K Cond a ks (S s - F s )M k ) t x (4) 

S=1 


where k = 0, 1,2. In this work, a c and K Cond are assumed 
to be unity. The polydisperse Schmidt numbers ( Sc k ) 
are given in Ref. [19]. In general, the coagulation and 
condensation/evaporation terms take different forms for 
continuum, slip, and free molecular regimes. One set of 
equations is used for all regimes by using a harmonic aver¬ 
age of the continuum and free molecular equations. 12 

A separate continuity equation is required for each gas 
phase species. Thus, for all species s, the saturation equa¬ 
tion takes the form 


dpC s 

dt 


+ v- 




1 

S Cs Re, 


V(ln C.) 


= Am, 


( 5 ) 


where C, is the gas phase mass fraction for species s. 
The RHS represents the sources in mass due to evapora¬ 
tion, condensation and chemical reaction. Coupling source 
terms for mass transfer between the bulk flow, individ¬ 
ual gas phase species and aerosol moments are determined 
from the relevant system parameters via 


Ns Ns 

Am = J2( Am s) = LX-auCS, - F s )M r ) (6) 

5=1 5=1 


where Am appears in the bulk flow continuity and Am, 
appears in each species equation. 

Flow, species and aerosol behavior are computed in 
separate modules using a finite difference discretization 


scheme and simultaneously marched in space based on the 
principles of the RNS methodology. 20 Multiple sweeps are 
required to converge the axial diffusion and axial pres¬ 
sure terms. Details of the discretization of the governing 
equations, implementation of boundary conditions, numer¬ 
ical coupling, matrix inversion routines and the local and 
global solution procedure are described in. 1718 

4. COMPUTATIONS OF GAS AND 
PARTICLE DYNAMICS IN MODEL 
CNT REACTOR 

In order to assess the impact of catalyst particle source 
size and concentration as well as wall conditions on cata¬ 
lyst particle dynamics in the CNT growth window, a series 
of computations have been carried out. Case 1 investigates 
overall effect of reactor conditions (temperature, veloc¬ 
ity and iron concentration) on the direction of the driv¬ 
ing force for particle evaporation/condensation. Cases 2 
and 3 additionally consider inlet and wall saturation con¬ 
ditions together with the effects of evaporation kinetics. 
Cases 4 through 6 investigate the influence of catalyst par¬ 
ticle concentration. Case 7 addresses the effect of C0 2 
production at the wall on the gas composition in the CNT 
growth window. Boundary conditions for each of the cases 
described below are listed in Tables I and II. Physical 
constants used in the calculations are listed in Table III. 
All cases were calculated in axisymmetric two dimen¬ 
sional space and included buoyancy with gravity pointing 
in the negative axial direction (i.e., flow from bottom to 
top). Ranges for source catalyst particle geometric number 
mean diameters (between 1.78 nm and 3.00 nm), concen¬ 
trations (between 1 x 10 14 #/m 3 and 1 x 10 18 #/m 3 ) and 
geometric standard deviation of the particle size distribu¬ 
tion (1.2) were estimated based on experimental observa¬ 
tions of collected samples 24 and particle size distribution 
measurements based on electrical mobility downstream of 
the reactor. 25 

4.1. Casel: Base Calculation of Gas Flow and 
Heat Transfer 

An initial single phase CFD computation was used to 
investigate the overall conditions in the reactor and, in 
particular, the direction of the driving force for particle 


Table I. Reactor conditions and gas phase boundary conditions for cases 1-6. 


Boundary 

condition 

Velocity (m/s) 

CO mole 

fraction 

N 2 mole 
fraction 

H 2 mole 
fraction 

Pressure (Pa) 

Temperature (K) 

HWG Inlet 

0.104793 

0.465 

0.465 

0.07 

Calculated 

298.15 

Coflow Inlet 

0.026947 

0.5 

0.5 

0.0 

Calculated 

298.15 

Reactor Wall 

0.0 

dC/dn = 0 

dC/dn = 0 

dC/dn = 0 

Calculated 

Set from 

experimental data 

Reactor Outlet 

dV/dn = 0 

dC/dn = 0 

dC/dn = 0 

dC/dn = 0 

1.01325 x 10 5 

dHt/dn = 0 

a dx/dn = 0 implies a 

zero gradient boundary conditions. 
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Table II. Catalyst particle boundary conditions for cases 2-6. 

Boundary 

condition 

Case 

Number 

concentration 

M 0 (#/m 3 ) 

Geometric 

mean diameter 

Dp (nm) 

Standard 

deviation 

of PSD 
<r g (N.D.) 

Fe vapor 
saturation 
ratio S (N.D) 
at reactor wall 

Catalyst 

i 

0.0 

0.0 

N.A 

1.0 

Particle 

2 

1 x 10 14 

3.0 

1.2 

1.0 

Source 

3 

1 x 10 14 

3.0 

1.2 

0.0 

Outlet 

4 

1 x 10 14 

1.78 

1.26 

0.0 


5 

1 x 10 16 

1.78 

1.26 

0.0 


6 

1 x 10 1S 

1.78 

1.26 

0.0 


7 

0.0 

0.0 

N.A 

0.0 

Coflow Outlet 

1-7 

0.0 

N.A 

N.A 

Calculated 

Reactor Outlet 

1-7 

dM 0 /dn = 0 

Calculated 

Calculated 

Calculated 


Table III. 

Properties of species used in calculations. 





Value 

Fe 

N, 

h 2 

C0 2 

CO 

Unit 

MW 

55.847 

28.013 

2.016 

44.010 

28.010 

g/mole 

y 

N.A 

1.40 

1.41 

1.28 

1.398 

ND 

Pr 

N.A 

0.71 

0.71 

0.71 

0.71 

ND 

Pr 

7040 

N.A 

N.A 

N.A 

N.A 

Kg/m 3 

S = A 

A = 2.41743 

N.A 

N.A 

N.A 

N.A 

J/mol K 

+ B T 

B =-3.1363 x 10“ 4 






+ C T 2 

C = 2.5893 x 10~ 9 






+ D T 3 

D = 2.7590 x lO" 13 






Cp = A 

A = 27.06 

A = 25.845 

A = 33.6667 

A = 42.388 

A = 25.6942 


+ B/T 2 

B = 0 

B = 47861 

B =-7307.16 

B = -889108 

B = 110861 


+ C/T 

C = 0 

C = 0.00 

C = -752.260 

C = 0.00 

C = 0.00 


+ D T 

D = 0 

D = 0.0089004 

D = -0.0109767 

D = 0.0150998 

D = 0.008292374 


+ E T 2 

E = 0 

E = -2.3480 x 10~ 6 

E = 1.2322 x lO' 5 

E = -2.9078 x 10~ 6 

E=-1.4769 x 10~ 6 


+ F T 3 

F= 0 

F = 2.1164 x lO” 10 

F = -4.7198 x 10- 9 

F = 0.00 

F = 0.00 


+ GT 4 

G = 0 

G = 0.00 

G= 8.1436 x 10- 13 

G = 0.00 

G = 0.00 


+ H T 5 

H = 0 

H = 0.00 

H =-5.2840 x 10~ 17 

H = 0.00 

H = 0.00 



(11.04482-18160.05715/T 
-1.85945 x 10 6 /T 2 






Ps 

10 +3.51941 x 10 S /T 3 ) 

N.A 

N.A 

N.A 

N.A 

Pa 

(T 

4.30 

3.798 

2.827 

3.941 

3.690 

m- 10 


3000 

71.4 

59.7 

195.2 

91.7 

K 


evaporation in the absence of the influence of catalyst 
particles. In general, since particle mass loadings were low 
for all the later cases investigated, these flow field and 
temperature calculations are also valid for cases 2-7. Iron 
vapor concentration and saturation ratios, however, change 
considerably under different catalyst particle source con¬ 
centration and wall saturation conditions. Figure 2 shows 
the conditions in the reactor under conditions in Table I, 
Case 1 where the walls are saturated with iron and parti¬ 
cle effects on the gas phase are neglected. Flow is from 
left to right with gravity pointing in the negative axial 
direction. Temperatures were found to be uniform and the 
entire flow was heated to within 5 K of the peak reactor 
temperature (Fig. 2(a)). The hot reactor section extended 
from approximately 0.45 m to 0.8 m. The flow in the 
reactor was found to be without recirculation up to and 
including the observed CNT growth window between 1167 


and 1179 K in the heating section (Fig. 2(b)) correspond¬ 
ing to a centerline location at approximately x = 0.39 
(0.14 meters from the HWG exit). Figure 2(c) shows the 
iron vapor concentration in the absence of iron catalyst 
particles and with iron saturated reactor walls. Iron vapor¬ 
ized from the walls was found to be relatively uniform 
across the reactor and reached a peak vapor mass frac¬ 
tion of approximately 1.1 x 1 0“ 7 near the hottest section 
and then decreasing as excess iron condensed on the 
reactor walls. By the reactor exit, the iron vapor con¬ 
centration dropped below 1 x 10 -9 . The effective satura¬ 
tion ratio for a given particle size including the Kelvin 
effect (taking into account the particle surface curvature) is 
defined as 


S* 



2S*v* " 
- r * k “ T f - 


( 7 ) 


J. Nanosci. Nanotechnol. 8, 1-17, 2008 


5 


RESEARCH ARTICLE 















RESEARCH ARTICLE 


CFD-Aerosol Modeling of the Effects of Wall Composition and Inlet Conditions on CNT Catalyst Particle Activity 


Brown et al. 



1406 

1332 

1258 

1184 

1111 

1037 

963 

889 

815 

741 

667 

583 

520 

446 

372 



0.372 

0.343 

0.314 

0.284 

0.255 

0.226 

0.196 

0.167 

0.138 

0.108 

0.079 

0.050 

0.021 

-0.009 

-0.038 


1.0E-07 

9.3E-08 

8.7E-08 

8.0E-08 

7.3E-08 

6.7E-08 

6.0E-08 

5.3E-08 

4.7E-08 

4.0E-08 

3.3E-08 

2.7E-08 

2.0E-08 

1.3E-08 

6.7E-09 



(d) 

0.04 

0.03 

0.02 

0.01 

I 0 

- 0.01 


9.00 

7.71 

6.43 

5.14 

3.86 

2.57 

1.29 

0.00 

-1.29 


-7.71 

-9.00 


Fig. 2. Gas temperature (a), velocity vectors shaded by axial velocity (b), iron vapor mass fraction (c), and effective saturation ratio for 3 nm iron 
particles (d), in the modeled reactor with no introduced pre-made particles at a wall peak temperature of 1480 K with coflow and particle source flow 
rates of 400 cm 3 /min. 


S* less than unity implies that particles will tend to evap¬ 
orate. Figure 2(d) shows an example of the calculated 
effective saturation ratio for 3 nm particles based on the 
expectation that evaporation would drive particles toward 
the observed diameter of 1.78 nm. Smaller particles have 
a higher effective saturation pressure and thus a stronger 
driving force for particle evaporation. Even for these rel¬ 
atively large catalyst particles, S* is less than 1 from 
approximately x = 0.41 to x — 0.85 m. It is apparent that, 
even in the presence of iron saturated walls there is a driv¬ 
ing force for particle evaporation and particles will tend 
to be reduced in size in the majority of the reactor. This 
calculation, however, does not take into account the kinet¬ 
ics of particle evaporation that would tend to reduce this 
driving force. Under the same flow rates and reactor tem¬ 
perature as in Case 1, additional cases were computed 
according to conditions in Table II in order to estimate 
these kinetic effects. Note that for clean (S — 0) reactor 
walls, there is no additional source of iron vapor and there 
always a driving force for catalyst particle evaporation. For 
brevity, results are not shown here. 

4.2. Cases 2-3: Effect of Wall Saturation Conditions 
Catalyst Particle Evaporation 

In Cases 2-3 the particle geometric number mean size, 
number concentration and standard deviation of the size 
distribution where kept constant while the wall concen¬ 
tration of iron vapor was varied in order to study the 
effects of wall saturation conditions on the resultant cat¬ 
alyst cluster dynamics. Calculations of the particle trans¬ 
port and dynamics were performed downstream from the 


source outlet plane at x — 0.25 until the reactor outlet. We 
assume a continuous change of particle size when, in fact, 
the evaporation process is discreet (atom by atom). As 
a lower bound, however, particle evaporation was halted 
when geometric mean diameter reached that of two iron 
atoms. Again, a source particle mean diameter of 3 nm 
was chosen. A source particle number size distribution 
of 1 x 10 14 #/m 3 was estimated based on experimentally 
observed concentrations at the outflow based on aerosol 
measurements. Note that the particle source concentra¬ 
tions are absolute values and are not normalized. A source 
particle geometric standard deviation of the PSD of 1.2 
was measured from High Resolution Transmission Elec¬ 
tron Microscopy images. 2 

Calculations showed that, at these conditions, particle 
concentrations were sufficiently low for inter-particle col¬ 
lisions and particle-to-gas coupling to be unimportant. 
Figures 3(a-c) shows the evolution of particle number 
concentration, mass fraction and geometric number mean 
size at the iron saturated wall condition. In Figure 3(a) 
we can see that the particle number concentration initially 
drops due to the reduction of gas density in the heating 
zone. The concentration then rapidly decreases at approx¬ 
imately x = 0.51 as clusters evaporate down to di-atoms 
and then rapidly diffuse to the reactor walls. In Figure 3(b) 
we can see the simultaneous effects of particle diffusion 
and evaporation. As the particles move downstream, there 
is initially diffusion in the radial direction which spreads 
the particle stream. At approximately x — 0.44 m, par¬ 
ticle begin to evaporate at an increasing rate due to the 
increased relative surface area and stronger Kelvin effect. 
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Fig. 3. Particle number concentration (a), mass fraction (b), and geometric number mean particle diameter (c), contours for 3 nm pre-made particles 
at a source concentration of 1 x 10 14 #/m 3 and a PSD geometric standard deviation of 1.2. Walls are saturated with iron. 


Figure 3(c) shows the evolution of particle diameter in the 
reactor. Here the simultaneous effects of preferential dif¬ 
fusion and evaporation can be seen. Due to their higher 
diffusivity, smaller particles diffuse more rapidly in the 
radial direction than larger particles. Thus the average par¬ 
ticle size is smaller near the reactor wall than at the reactor 


difference in the location of particle evaporation) indicat¬ 
ing wall saturation conditions do not play a significant role 
in particle evaporation dynamics and thus are not impor¬ 
tant in explaining either the initiation or the termination of 
CNT synthesis in the observed CNT growth window. 


core. Rapid particle evaporation occurs between x = 0.41 
and 0.51 m. 

Overall, particles were found to be stable up to approx¬ 
imately 0.18 m from the source outlet (x = 0.43 m) and 
did not evaporate until approximately 0.26 m downstream 
(x = 0.51 m). All particles survived with no apprecia¬ 
ble reduction in size up to the experimentally measured 
CNT growth window 0.14 m from the source outlet 
(x = 0.39 m). Lowering the iron vapor concentration at 
the reactor walls (simulating a clean ceramic reactor tube) 
produced virtually identical results (less than a 1 mm 


4.3. Cases 4-6: Effect of Inlet Catalyst Particle 
Concentration on Particle Dynamics 

In order to determine the effect of catalyst particle concen¬ 
tration on CNT growth dynamics, Cases 4-7 were com¬ 
puted by keeping the wall saturation conditions as well 
as particle geometric mean diameter and standard distri¬ 
bution of the particle size distribution (1.78 nm and 1.26 
respectively) constant and varying the particle source con¬ 
centration from 1 x 10 14 to 1 x 10 18 #/m 3 . Figure 4 to 
Figure 7 show the dramatic effect of increasing the inlet 
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Fig. 4. Particle number concentration along the centerline (a), of the reactor at a wall peak temperature of 1480 K for source particle concentrations 
of 1 x 10 14 , 1 x 10 16 and 1 x 10 ls #/m 3 and particle number concentration contours for each inlet concentration (b, c, and d). 
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Fig. 6. Particle mass fraction along the centerline (a), of the modeled reactor at a wall peak temperature of 1480 K for source particle concentrations 
of 1 x 10 14 , 1 x 10 16 and 1 x 10 18 #/m 3 and particle mass fraction contours for each inlet concentration (b, c, and d). 
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Fig. 7. Standard deviation of the PSD along the centerline (a), of the modeled reactor at a wall set temperature of 1200 C for source particle 
concentrations of 1 x 10 14 , 1 x 10 16 and 1 x 10 18 #/m 3 and standard deviation of the PSD contours for each inlet concentration (b, c, and d). 
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particle number concentration on the catalyst particle size 
and size distribution in the reactor. In all cases, particles 
are stable up to and including the observed CNT growth 
window at x = 0.39 m, however, higher particle concen¬ 
trations result in significant particle collisions and cause 
the particle number concentration to decrease (Fig. 4). 
Higher initial concentrations lead to more rapid reduc¬ 
tion rates (Fig. 4(a)). At the highest source concentration 
(1 x 10 18 #/m 3 ), the catalyst particle number concentration 
decreases by more than 3 orders of magnitude by the time 
particles are in the CNT growth window and are actually 
lower in concentration than if they had been introduced 
at a number concentration two orders of magnitude lower 
(1 x 10 16 #/m 3 ). The evolution of particle size is shown 
in Figure 5. Interparticle collisions and coagulation causes 
particles to grow and, at higher concentrations, the parti¬ 
cles grow more rapidly (Fig. 5(a)). At the highest calcu¬ 
lated source concentration, the particles are found to grow 
by almost a factor of 10 (from 1.78 to 15.0 nm) by the 
time they reach the CNT growth window. Such particles 
are generally too large to synthesize SWCNTs. As can 
be seen in Figure 5(b), 1.78 nm particles do not readily 
collide and coagulate into larger particles at a source con¬ 
centration of 1 x 10 14 #/m 3 , but do diffuse more rapidly 
than the 3 nm particles of Cases 2 and 3 leading to a more 
even radial distribution in the reactor. They are found to 
evaporate approximately 0.05 m earlier than the 3 nm par¬ 
ticles. Figure 5(c) shows that 1.78 nm particles introduced 
at a concentration of 1 x 10 16 #/m 3 grow to approximately 
3 nm near the centerline before evaporation begins and to 
approximately 2 nm near the wall. These particles evap¬ 
orate at approximately x = 0.51 m. At the highest source 
concentration, particles grow large enough not to evapo¬ 
rate and are able to exit the reactor intact (Fig. 5(d)). This 
is due to a reduced Kelvin effect, a reduced ratio of sur¬ 
face area to particle mass and an increase in the iron vapor 
mass fraction due to evaporation from particles. Addition¬ 
ally, these particles have very low diffusivity and remain 
close to the reactor centerline in a well defined column. 
The combined effects of particle collisions on diffusion 
and evaporation are seen more clearly in Figure 6 showing 
that higher source concentrations lead to larger, less dif¬ 
fusive and more slowly evaporating particles. Not only do 
high concentrations of source particles increase the mean 
particle size, leading to a large fraction of “dead” catalyst 
particles that are too large to synthesize CNTs, but it also 
leads to a significant widening of the particle size distri¬ 
bution as can be seen in Figure 7, resulting in less con¬ 
trol over CNT diameter and, consequently, chirality. At the 
highest source concentration (lx 10 18 #/m 3 ), the geometric 
standard deviation rises to the theoretically limiting value 
of 1.355 approximately 0.05 m from their point of intro¬ 
duction. At the intermediate source concentration (1 x 10 16 
#/m 3 ), the geometric standard deviation approaches but 
does not reach this limiting value. This explains why the 


highest calculated source concentration (1 x 10 18 #/m 3 ) 
leads to a lower number concentration in the growth win¬ 
dow: wider PSDs have increased collision efficiencies. 

4.4. Cases 7: Prediction of C0 2 Production Due to 
Catalytic Wall Reactions 

As is apparent from the previous calculations, iron cat¬ 
alyst particles are not expected to evaporate before the 
growth window under even the most severe reactor con¬ 
ditions and the lowest particle loadings. Consequently, the 
effect of catalytic wall reactions was investigated as an 
explanation of the experimental observations of CNT syn¬ 
thesis. C0 2 in the gas phase has been shown to promote 
the synthesis of CNTs. 26 The importance of having iron 
on the reactor walls can be explained by its catalytic activ¬ 
ity in CO disproportionation and hydrogenation reactions 
that lead to the release of gaseous products, such as car¬ 
bon dioxide (C0 2 ) and water (H 2 0). It was shown that 
those gaseous components formed upstream on the reactor 
walls play an important role downstream in the region of 
the CNT growth. This hypothesis is supported by the work 
of Hata et al. 27 and Nasibulin et al., 14 who have recently 
demonstrated the importance of water vapor and C0 2 in 
providing conditions for efficient growth of CNTs. Thus, 
an additional series of calculations was carried out in order 
to determine the correlation of catalytically produced C0 2 
to the location of the growth window. CO/C0 2 conversion 
data 13 (Fig. 8) showed a kinetically limited region up to 
approximately 1170 K according to the reaction 

CO(g) + CO(g)^C0 2 (g) + C(s) (8) 

and a thermodynamically limited region above approxi¬ 
mately 1,170 K. From this data, we have calculated an acti¬ 
vation energy (Ea) of 63 kJ/mole in the kinetic region and 
this reaction was incorporated as CO and C0 2 boundary 

1171 (K) Ln(X C o 2 ) = 7.16091-7581.28855/T 



1/T (1/K) 


Fig. 8. Conversion of CO to C0 2 versus temperature showing kineti¬ 
cally limited and thermodynamically limited regions. 
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Fig. 9. C0 2 mass fraction contours in the modeled reactor showing catalytic production of C0 2 at the wall upstream of the observed CNT growth 

window. 


conditions in the calculations such that, below the ther¬ 
modynamic equilibrium state, CO conversion was limited 
kinetically according to reaction (8) with the rate of pro¬ 
duction of C0 2 is given by 


^ Ox ) 2 

d t 


= A f Exp 

S 


Ea 

ItT* 


C, 


co 


( 9 ) 


Note that the reaction rate is then per unit area of the 
reactor wall. In the thermodynamics equilibrium region, 
the measured equilibrium concentration was used for the 
CO and C0 2 boundary conditions. Additionally we did 
not allow any catalytic decomposition reactions (kinetic 
or equilibrium) at walls where iron particles would be 
expected to be unstable. This lower limit was estimated 
to be 1200 K based on the earlier particle evaporation 
calculations. Experimental observations did not allow the 
accurate determination of the forward reaction rate A*, 
however it was experimentally observed that, at the outlet 
of the reactor, the concentration of C0 2 was approximately 
150 ppm as measured by Fourier Transform Infrared Spec¬ 
troscopy (FTIR) when iron was present on the reactor 
walls. CFD calculations showed that assuming A* = 2.0 x 
10“ 3 mole _1 m _2 s _1 gave similar outlet concentrations 
(Fig. 9). We did not consider the yield of H 2 0 since FTIR 
measured values were an order of magnitude lower than 
for C0 2 . 28 Figure 9 shows the mole fraction of C0 2 in the 
reactor clearly indicating that the location of the growth 
window approximately 0.18 m from the catalyst source 
outlet (x = 0.43) corresponds to a significant product con¬ 
centration (approximately 80 ppm) and indicates a strong 
correlation between CNT formation and growth and the 
local availability of C0 2 . 

As for the role of C0 2 in the reactor, several theories 
have been introduced. Adier et al. 28 reported C0 2 to be 
important to prevent the conversion of iron (Fe) to cemen- 
tite (Fe 3 C) which is likely to poison CNT catalyst particles. 
Nasibulin et al. 28 additionally reported the essential role of 


C0 2 for the production of CNTs and, specifically, proposed 
it as an oxidizing agent that etches the surface of catalyst 
particles and CNTs and thus removes amorphous carbon 
which can block the introduction of additional carbon or 
poison the surface of catalysts and CNTs. 

5. DISCUSSION 

Calculations of the flow, heat transfer and gas and particle 
phase behavior of iron in the CNT reactor have revealed 
the insensitivity of the resulting catalyst particle size to 
wall Fe saturation conditions and the strong sensitivity 
of the resulting catalyst particle size and polydispersity 
to their initial size and concentration. In all cases and 
under the most extreme temperature conditions, even the 
lowest particle source concentrations (1 x 10 14 #/m 3 ) and 
smallest sizes (1.78 nm), catalyst clusters were found to 
be stable up to the observed CNT growth window. How¬ 
ever, high concentrations resulted in significant particle 
growth by collision, effectively deactivating a large frac¬ 
tion of particles and widening their distribution. The result¬ 
ing wider particle size distribution and increased mean 
diameter actually leads to an overall reduction in catalyst 
particle concentration in the growth window at the highest 
source concentrations calculated (1 x 10 18 #/m 3 ) as com¬ 
pared to the intermediate concentration (1 x 10 16 #/m 3 ). 
Consequently, large catalyst particle concentrations can, 
paradoxically, lead to lower yields of CNTs of lower 
homogeneity. Inflow concentrations should thus be kept 
below 1 x 10 16 #/m 3 to maintain favorable properties for 
controlled CNT production. The composition of the reac¬ 
tor wall was found to be important, not due to the sup¬ 
ply of iron vapor, as was originally believed, but, instead, 
due to its ability to supply C0 2 (through catalytic wall 
decomposition reactions) which is understood to enhance 
CNT nucleation and growth. Consequently, the likely role 
of the conditioning of the reactor wall was found to be 
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chemical (as a source of CNT precursor and/or promoter 
such as C0 2 ) rather than physical (as a source of sufficient 
Fe vapor to prevent evaporation). As the catalyst diame¬ 
ter in the active growth window is the critical factor in 
determining the diameter of synthesized CNTs, it is, there¬ 
fore, vital to simultaneously control the concentration and 
diameter of catalyst particles and the concentration of C0 2 
in important regions of the reactor in addition to the CO 
concentration and gas temperature. 

It is useful to interpret these results in terms of the avail¬ 
able literature for floating catalyst CNT production. For 
instance, in the HiPCO process, 4 catalyst particles are gen¬ 
erated in situ via the thermal decomposition of Fe(CO) 5 
vapor under high pressure conditions in a quartz reactor. It 
is known that a large fraction of catalyst particles remain 
inactive, that the distributions of catalyst particle size and 
produced CNT chirality are wide and that there is a sig¬ 
nificant amount of amorphous carbon remaining on the 
produced CNTs. This can be explained in terms of the cat¬ 
alyst particle production and collision dynamics and the 
availability of reagents such as C0 2 . In terms of the parti¬ 
cle dynamics, in the HiPCO process, iron catalyst particles 
are produced continually in a wide reaction zone where 
the precursor gas decomposes to release supersaturated Fe 
vapor. This results in a large reaction zone where compe¬ 
tition between nucleation and collision and high particle 
number concentrations tend to broaden their size distribu¬ 
tion. Moreover, our results indicate that, since the HiPCO 
reactor has significant exposed surface area at moderate 
temperatures, one can expect the release of C0 2 from the 
reactor walls in the same manner as proposed in our reac¬ 
tor. In fact, as was reported in (Ref. [29]), a high con¬ 
centration of C0 2 (on the order of 10,000 ppm) is known 
to exist at the HiPCO reactor outlet under CNT synthesis 
conditions. Though they have not reported C0 2 concen¬ 
trations in the absence of Fe(CO) 5 precursor, the presence 
of reactor walls and turbulent mixing (which leads to the 
deposition of iron on the walls) provides additional sources 
of C0 2 than from conversion on floating catalyst parti¬ 
cles alone. Additionally, at high pressures the equilibrium 
of CO disproportionation reaction (1) is shifted towards 
the formation of C0 2 and C. 30 Thus, the high CO pres¬ 
sure effect is not only in the enhancement of the reaction 
rate due to an increase of the amount of reagent, but also 
in widening of the CNT formation temperature range and 
enlarging the growth window so that CNT growth occurs 
even up to 1500 K at 30 bars. 30 Thus, the high concen¬ 
tration of inactive catalyst particles having a wide size 
distribution in the HiPCO process can be explained as a 
result of the insufficient control of particle size due to 
their high concentration and continuous nucleation. Addi¬ 
tionally, C0 2 concentration (which can inhibit coagulation 
and growth of active catalyst particles by initiating CNT 
growth) is unknown in the widened growth window. 

Zhao et al. 31 noted the importance of reactor conditions 
on the efficiency of catalyst use for CNT production from 


sulfur contaminated ferrocene vapor. They reported that, 
for their experimental setup, reducing the diameter of the 
growth tube resulted in more efficient use of catalysts. 
This can be explained in terms of our current results in 
that smaller reactor tubes led to higher flow rates, and 
shorter residence times in the reaction zone. Shorter resi¬ 
dence times further reduced the time for precursor decom¬ 
position and particle-particle interaction, thus maintaining 
a smaller and possibly narrower particle size distribution 
upon CNT nucleation, thus increasing yield efficiency. 

Zheng et al. 32 produced single-walled CNTs by a CVD 
method wherein they introduced CO gas pretreated in 
a 773 K furnace. They noted that the pretreatment furnace 
“was crucial to obtain clean single-walled CNT sample.” It 
was explained that the pretreatment was needed to remove 
Fe(CO) 5 impurities in commercial CO. In the framework 
of our computations, however, we feel their result indicate 
that the pretreatment stage was necessary both to decom¬ 
pose Fe(CO) 5 on the pretreatment reactor walls and to 
release C0 2 due to the CO disproportionation reaction on 
the deposited iron. Thus, their pretreatment reactor 32 and 
our heated walls upstream of the growth window serve to 
provide CNT growth promoters (in this case C0 2 ) known 
to aid the synthesis of CNTs. 

6. CONCLUSIONS 

A CFD based polydisperse aerosol population balance 
model suitable for predicting the behavior of iron cata¬ 
lyst particles and surface reactions is used to understand 
the effects of catalyst particle size and concentration, sur¬ 
face CO disproportionation chemistry, heat transfer and 
wall saturation levels on catalyst particle size and synthesis 
conditions in floating catalyst synthesis methods. In this 
work, the dynamics of preformed catalyst particles intro¬ 
duced near the centerline of a heated tubular reactor into 
conditions leading to CNT nucleation and growth was sim¬ 
ulated. It is found that introduced catalyst particles with 
diameters suitable for catalyzing SWCNTs do not evapo¬ 
rate until CNT nucleation under all studied reactor condi¬ 
tions. In particular, it is found that the level of wall iron 
vapor concentration does not significantly affect catalyst 
particle evaporation. However, catalyst particle concentra¬ 
tions significantly above 1 x 10 14 #/m 3 lead to poor control 
over particle size and polydispersity at the CNT nucleation 
front which, in turn, is expected to lead to poor control 
over produced CNTs. The location of the experimentally 
observed growth window is shown to be associated with 
the availability of C0 2 produced by catalytic wall reactions 
and not with the stability of catalyst particles. Together, 
these results demonstrate the necessity of simultaneously 
controlling CO, C0 2 and catalyst particle size and concen¬ 
tration as well as gas temperatures and residence times in 
critical regions of the reactor in order to synthesize well 
controlled CNTs. 
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Ath volume moment of the distribution 


9. 


M t = f v *p n *p ( v * p > t*) dv *r = M o v f ex p(Jk 2 1 h 2 or 

/ Ns 

Mixture molecular weight I y ~^p s C s + p n C, 


PsC s 

h MW , 


Particle number concentration 
Number of gas phase species 
Pressure 

Prandtl Number (Cp/i / /k) 

Heat flux 

Universal gas constant 

Laminar Reynolds number (/kL 0C L 0C //x, 0 J 

Local Reynolds number (p* U^L^/p,*) 

Surface area 

Saturation ratio 


(ND) 

(ND) 

(ND) 

(ND) 

(m/s) 

(m 3 /s) 

(m 2 /s) 

(m 3 /s) 

(ND) 

(ND) 


Cp 

Mixture specific heat at constant pressure I y^ C,Cp, + C„ Cp„ 1 

(J/kg-K) 

d 

\S=1 ' 

Diameter 

(m) 

D 

Diffusion Coefficient (k B Tfr* / m .*) 

(m 2 /s) 

Ect 

Activation Energy 

(J/mole) 


(ND) 

(ND, m/s 2 ) 
(m 3 /m 3 /s) 
(m 2 /s 2 ) 
(mole/m 2 /s) 
(kg m 2 /s 2 /K) 
(ND) 
(ND) 
(ND) 

(m) 

(kg) 

(ND) 

(m 3 *7m 3 ) 


(kg/mol) 

(#/m 3 ) 

(ND) 

(kg/m/s 2 ) 

(ND) 

(kg/s 3 ) 

(kg m 2 /s 2 /Kmol) 
(ND) 

(ND) 

(m 2 ) 

(ND) 
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Sc 


T 

u, v, w 

u, y, w 
u 

v 


X, y,z 

<*k 


P 

y 

y k 


A 

e* 

K 

A 

</> 

M 

P 

cr 

2 

T 


Schmidt number ( Sc = 


P-] 


Pf 18 m; 


P/Dp pfk B T*dfC 


Sc 0 = l/(A 3 v“ 1/3 exp(0.5 In 2 cr) + A 4 v“ 2/3 exp(2 In 2 cr)) 

Sc { = l/(A 3 v“ 1/3 exp(—2.5 In 2 cr) + A 4 v“ 2/3 exp(—4In 2 cr )) 

Sc 2 = l/(A 3 v“ 1/3 exp(—5.5 In 2 cr) + A 4 v“ 2/3 exp(—10In 2 cr)) 
Temperature 

Cartesian components of velocity 

Contravariant components of velocity 

Velocity vector 

/ M, 2 

volume v„ = —^^ 
s M n L 5 M ?' 5 


Cartesian coordinates in physical space 
Ath condensation coefficient 


= 0, a u = 


c 




= 


c 


z,+<£\U' 


a u = V /3 exp(-41n 2 cr) a["‘ = B ls v g 1/3 exp( -- In 2 cr 


11 


a 2 s = B 3 s v s ’ exp(-lOln-cr) a£ = B ls v g 1/3 exp( - — In cr 

Collision frequency function 
Mixture ratio of specific heats (1 — R/MW Cp) 
kth coagulation coefficient 


To = 


c fwi 

To To 


To + To 


M 0 X L X 


Tt =0, y 2 = 


Hy? 

y^y f 2 


M 2oo L oa 


To = ~B a 


1 


1 + exp(ln’ cr) + B 5 Kn exp I - In" cr ) (1 +exp(21n cr)) 


To 


= -B,bv r y 6 


25 


8 


exp I — In’ cr 1 + 2 exp I - In" cr ) + exp I - In cr 


8 


1 


8 


2B a 


y 2 = — exp(-271n- cr) 


(1 + exp(ln 2 s)) + B s Kn exp( — - In 2 cr ) (1 + exp(—21n 2 cr)) 


T2 


fm 


2 B 2 b 2 

/ 11/6 

Vg 


exp(—36 In 2 cr) 


exp 


85 

¥ 


In 2 cr )+ 2 exp I — In 2 cr 


8 


- exp 


109 


In" cr 


Change due to aerosol formation and growth 
Lennard-Jones Temperature 
Thermal conductivity 


. P'f TTm* 

Mean free path of gas mixture A* = — 1 


2 \ + d„ 


Pf 


2k b T* 


Knudsen correction , 

v d p + 5.33(A 2 /c/ ;i ) + 3.42A 

Viscosity 

Mixture density Ps + Pa 

Standard deviation of aerosol size distribution 

In 2 cr = - In 
9 

Surface tension 


M 0 M 2 

Ml 


+ hr cr^ ), Lennard-Jones Diameter 

py?c c 


Viscous stress tensor, characteristic time r = 


18 iL* 


Transformed coordinates in computational space 


(ND) 


(K) 

(m/s) 

(ND) 

(ND) 

(m 3 ) 

(m) 

(ND) 


(#/m 3 /s) 

(ND) 

(ND) 


(ND) 

(K) 

(kg/s 3 /K) 

(m) 

(ND) 

(kg/m/s) 

(kg/m 3 ) 

(ND, A) 


(ND, N/m) 
(kg/m/s 2 , s) 
(ND) 
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Subscripts 


References and Notes 


a aerosol phase 

j L 17 , c 

Con due to condensation 

C convective 

D diffusional 

Eq equilibrium 

/ bulk fluid 

g geometric mean 

k Ath aerosol moment 

/ laminar 

L local 

p particle 

s gas phase species s 

t total 

T thermophoretic 

v viscous 

w wall value 

2 surface tension 

7), £ differentiated with respect to 17 , £ 

00 initial, free stream or characteristic value 


Superscripts 

* dimensional quantity 

i x, y, z 

c continuum 

fm free molecular 

k to the power k 


Non-Dimensionalization 


uj 




p = 


Pi 

Poo 


p = 



t = 


LJUJ tx U x 


M k = 


Ml 

M k J 



4 

- 77 r 

3 poo 
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